Computer Methods in Biomechanics and Biomedical Engineering

Acute subdural hematoma (ASDH) is a prevalent injury with high mortality and morbidity, often resulting from bridging vein (BV) disruption secondary to cortical relative motion. As a thin membrane enveloping the brain surface and anchoring BVs, the pia mater is hypothesized to play a critical mechanical role in cortical response and hence ASDH pathogenesis. Finite element (FE) head models are valuable tools to predict ASDH occurrence during impacts. However, the pia mater is often represented as an elastic material in existing FE head models, despite experimental evidence reporting its nonlinear mechanical behavior. In this study, both linear (Young's modulus of 11.5 MPa) and nonlinear (the stress-strain curve derived from pial tension tests) material models of the pia mater were implemented in one FE head model. The models were subjected to three experimental impact loadings, one of which was known to cause ASDH and two of which were not. Results demonstrated that, across all simulated impacts, the model with nonlinear pia mater properties predicted larger cortical displacements and BV responses than the linear model. For the impact with known ASDH occurrence, the predicted BV strain was 0.17 for the nonlinear model and 0.094 for the linear model, with only the former approaching the reported rupture strain range of the BV-superior sagittal sinus complex (0.29 {+/-} 0.13). These findings verified the mechanical importance of the pia mater in cortical responses and hence the prediction of ASDH, suggesting that conventional linear pia modeling might over-constrain cortical motion, leading to underestimation of BV strain and ASDH risk. The current study supported the adoption of experimentally derived nonlinear pia mater properties in FE head models to improve the reliability of ASDH prediction.

Trees are routinely pruned to mitigate the risk of wind damage, but there are few studies examining changes in wind loads after pruning, especially for large conifers. In this study, ten Colorado spruces (Picea pungens) were monitored before and after a series of pruning treatments. Trees were pruned to raise or thin crowns over a range of severities between 0% and 40%. Wind-induced bending moments were measured using two calibrated displacement probes installed orthogonally on the lower stem of each tree. Using a hierarchical Bayesian model, the relationship between maximum wind speeds and bending moments was quantified, consistent with theoretical and empirical expectations, as a non-linear power law. Random intercepts for model coefficients were used to account for individual variability in aerodynamic behavior among experimental trees, and predictions were made using the median response marginalized over the observed trees. The modeled relationship between wind speeds and bending moments was physically reasonable and like existing measurements with scaling exponents below two. Despite considerable variation among experimental trees, the aerodynamic behavior of trees, as indicated by model coefficients, was not clearly altered by pruning treatments, and, correspondingly, model predictions of bending moments over the range of observed wind speeds remained similar for all pruning treatments. Ultimately, the study yielded weak evidence for a change in bending moments following conventional pruning treatments for Colorado spruce, and the practical value of pruning to mitigate risk appeared limited for the studied conditions. Highlights- Wind loads were monitored on large Colorado spruce after crown raising and thinning - A hierarchical Bayesian model quantified wind speed and bending moment power laws - Negligible change in bending moments was found for all pruning types and severities - Conventional pruning methods may not mitigate risk for Colorado spruce

This study examined how tempo, dynamics, and string influence upper-extremity physical exposure in professional violinists and how exposure variability is distributed among musical characteristics, between-subject differences, and residual variability. Twelve violinists performed seven standardized scales while bilateral upper-arm and wrist kinematics and shoulder and forearm muscle activity were recorded. Linear mixed-effects models showed that faster tempo increased right upper-arm velocity and bilateral forearm activity while reducing right upper-arm and wrist ranges of motion. Louder dynamics increased bilateral forearm and right trapezius activity and right-wrist ranges of motion. Higher-posture strings increased right upper-arm elevation and right shoulder muscle activity. Variance analysis identified exposures predominantly related to musical characteristics, jointly related to musical characteristics and between-subject differences, predominantly related to between-subject differences, or mainly unexplained. These findings support future exposure prediction from musical characteristics and targeted prevention through repertoire-based workload management, structured recovery, and individualized technique-focused strategies.

The control of uterine activity during pregnancy is a complex process that involves regulating myometrial excitability across multiple scales. While numerous studies have investigated various regulatory mechanisms and established the contributions of ion channels and gap junctions, how these mechanisms interact to produce observed changes in uterine activity remains poorly understood. Pivotal to these efforts are computational models that effectively capture gestational changes in excitability across scales. In this study, we propose a multiscale computational modeling framework that can reproduce measured activity at the cellular and tissue scales at a given gestational stage. At the cellular level, we identify key ion currents underlying the observed electrophysiological properties based on a literature review of their regulation and a sensitivity analysis of the Tong 2011 uterine smooth muscle cell activation model. The conductances of these ion currents are then fit to reproduce characteristic resting membrane potentials and burst properties using Bayesian optimization. To extend to the tissue level, we employ an anisotropic monodomain model, parameterized by the resistivity of late pregnancy uterine muscle, to investigate electrical propagation in a two-dimensional section of uterine tissue. We then apply the multiscale model to study myometrial activation in late pregnancy and elucidate the contributions of ion channel and gap junction regulation in transitioning the uterus from a quiescent state to labor. Our resulting model successfully reproduces measured electrophysiological properties at the cellular level and characteristic single-spike and burst-propagation patterns at the tissue level across the three late-pregnant time points analyzed (days 16/17, 18/19, and 20/21) in a murine model. Furthermore, our results suggest that the regulation of the conductances of the voltage-dependent potassium current (IK1), L-type calcium current (ICaL), and sodium current (INa) is most important in determining preterm uterine excitability. The framework established here will promote the development of more gestationally relevant models to better understand labor progression and the factors involved in dysfunctional labor.

Muscle force sharing during locomotion is influenced by the mechanical demands of movement and the contractile properties of synergistic muscles. In cats, plantarflexor muscles exhibit distinct functional specialization, with the slow-fibred soleus maintaining relatively constant force across conditions while faster muscles such as the plantaris and gastrocnemius increase force production with increasing locomotor demand. However, it remains unclear whether similar force-sharing patterns occur in larger animals with different musculoskeletal designs. Therefore, the purpose of this study was to examine force sharing between the superficial digital flexor (SDF) and medial gastrocnemius (MG) muscles during treadmill locomotion in sheep. Tendon buckle force transducers were surgically implanted on the SDF and MG tendons of seven sheep, and in vivo muscle forces were recorded during walking and trotting across different speeds and inclines. Both muscles increased force with increasing speed and incline; however, speed had a substantially greater effect than incline. The SDF consistently produced greater absolute force than the MG across all conditions, whereas the MG exhibited slightly larger relative increases in force with increasing speed. Time to peak force decreased with increasing speed in both muscles, although the SDF reached peak force later in stance than the MG across conditions. In contrast to the distinct specialization observed in cats, neither muscle displayed a relatively condition-independent, soleus-like force contribution. These findings suggest that force sharing in sheep is more distributed across synergistic muscles and may reflect the influence of musculoskeletal design, tendon compliance, and mixed fibre-type composition on muscle function in larger species.

Total knee replacement restores mobility in patients with advanced osteoarthritis, yet many individuals still experience limited ability to perform high-flexion tasks such as squatting. Current preoperative planning relies on static imaging and cannot predict how different implant alignment choices will affect postoperative dynamic function. This study developed a predictive simulation framework that uses bi-level inverse optimal control to link preoperative implant alignment directly to expected postoperative squat kinematics. Subject-specific musculoskeletal models were constructed for six total knee replacement patients using experimental squat data. Bi-level inverse optimal control was applied to identify both individualised and group-level cost functions. The individualised setting provided subject-specific accuracy, while the group-level setting derived a single group-level cost function as an initial step toward preoperative use without requiring postoperative motion data. The individualised setting reproduced experimental trajectories with low errors across all joints (mean apex difference 1.53{degrees}, root-mean-square error 5.15{degrees}, normalised root-mean-square error 11.15%, Pearson correlation 0.96). The group-level setting yielded higher but acceptable errors (mean apex difference 5.70{degrees}, root-mean-square error 6.75{degrees}, normalised root-mean-square error 17.53%, Pearson correlation 0.95) while preserving the general pattern and phasing of the motion. Squat depth emerged naturally from the optimisation rather than being prescribed. This framework may provide a basis for future quantitative tools to evaluate how implant alignment choices influence postoperative squat performance, potentially improving functional outcomes in total knee replacement. These results suggest that the proposed IOC framework can reproduce key features of post-TKR squat kinematics, but further out-of-sample validation is required before it can be used for preoperative prediction or translated into tools aimed at improving functional outcomes in total knee replacement.

Computational Hill-type muscle models are widely used to simulate muscle force production because of their efficiency and physiological interpretability. However, their formulation relies on limiting assumptions, including debated multiscale simplifications, a simplified excitation-activation dynamics and an inability to capture slow and fast fibres. Moreover, existing Hill-type models remain insufficiently validated across physiological scales, fibre types, and contraction modes. We addressed these limitations by developing a multiscale fibre-type specific Hill-type neuromuscular actuator with mechanistic excitation-activation dynamics and systematically validated it against comprehensive experimental benchmarks. The model built upon a previously proposed motoneuron-driven actuator incorporating calcium-kinetics-based activation dynamics. The excitation-activation formulation was further refined to strengthen its physiological basis, while the contraction dynamics was extended by including an activation- and length-dependent force-velocity relationship, elastic tendon, passive elastic element, and the fibre-type-specific effects of yielding and sag. Validation was performed against four benchmark datasets spanning motor-unit and whole-muscle scales, including slow and fast fibres under both isometric and dynamic conditions. Experimental force traces were obtained from six muscles of rats and cats using a broad range of stimulation frequencies, muscle lengths, and imposed length changes, combining previous literature datasets with experiments performed ad hoc for this study. Overall, the model reproduced forces across all benchmark conditions, with mean absolute errors typically below 15% of the maximum isometric force, although larger errors were observed in specific submaximal and dynamic trials. The inclusion of physiologically based excitation-activation dynamics, together with yielding and sag, improved model performance under submaximal activation conditions. This study presents the first systematic validation of a single multiscale Hill-type neuromuscular actuator against comprehensive experimental motor unit and muscle force data, providing a benchmark framework for the development and assessment of future models. Author summarySkeletal muscles generate force through a complex sequence of events that links neural signals to muscle contraction. Because direct measurements are difficult to obtain, researchers often rely on computer models to investigate neuromuscular function and estimate muscle forces. However, most modelling approaches rely on simplifying assumptions about how force is generated across different biological scales, how muscles are activated, and how slow and fast muscle fibres behave. Moreover, they have not been validated against comprehensive experimental data. As a result, it remains unclear how accurately these models can reproduce muscle force across different physiological conditions. In this study, we established the first comprehensive set of experimental benchmarks spanning both motor-unit and whole-muscle scales, including slow and fast muscles under isometric and dynamic conditions. We used these benchmarks to validate a newly developed multiscale muscle model that explicitly represents the physiological pathway from neural stimulation to force production. The model incorporates experimentally based descriptions of calcium dynamics, activation, tendon elasticity, and fibre-type-specific contractile properties. We then compared simulated and experimental force responses across a wide range of stimulation frequencies, muscle lengths, and length-change conditions.

The osteocyte network in bone is believed to play an important role for how bone tissues sense and respond to mechanical stimulation. Yet, bone adaptation to mechanical loads is often conceptualised as a simple response to mechanical stimuli, such as Wolffs law, which is based on mechanical variables only and takes no account of the cellular basis of mechanosensation. Wolffs law presumes the existence of a reference mechanical stimulus, the mechanical setpoint, above which bone is consolidated, and under which bone is removed. In this paper, we develop a theory of bone tissue sensing and adaptation based on osteocytes to provide new understanding of the role played by osteocyte signals in mechanical adaptation. In this theory, the mechanical setpoint of Frosts mechanostat is explicitly embodied as osteocyte properties involved in mechanotransduction. The mechanical setpoint is allowed to adapt due to the replacement of osteocytes during remodelling, making the setpoint space and time dependent. We propose a mathematical model to implement this new theory of bone adapation and present numerical simulations of this model to explore how mechanobiological response curves (effective Wolffs laws) are modulated by setpoint adaptation during remodelling. By accounting for varying osteocyte populations within bone tissue, we explore bone adaptation under osteocyte disruptions, which is particularly relevant to age-related bone loss. Our model suggests that biological disruptions of remodelling balance cannot always be compensated by mechanical feedback, and that setpoint adaptation during remodelling may have significant observable consequences, such as hysteresis in bone response signatures that resemble lazy zones.

Objective: To identify brain injury metrics suitable for supporting sports head injury assessment by evaluating their association with brain tissue strain and consistency across sports. Methods: Head kinematics from 3,139 impacts in boxing, mixed martial arts, and rugby matches were recorded using instrumented mouthguards and used to calculate nine brain injury metrics. Impacts were simulated using an anatomically detailed finite element brain model to estimate peak 95th-percentile maximum principal strain (MPS) in the brain and brainstem, a measure of tissue deformation associated with long-term pathology. Sport-specific ordinary least squares models estimated xE, the metric value equivalent to a reference MPS of 0.21. Metric-MPS correlations and xE uncertainties were quantified using 5000 bootstrap resamples. Cross-sport consistency was assessed using the coefficient of variation (CV) of sport-specific median xE values, and uncertainty using the normalised confidence interval size (NCIS). Results: XGB, an extreme gradient boosting strain-prediction model, showed the strongest and most consistent correlations with whole-brain (r=0.924-0.974) and brainstem MPS (r=0.887-0.954) across all sports. PRV, BrIC and UBrIC also correlated strongly with whole-brain (r=0.724-0.930) and brainstem MPS (r=0.739-0.900), whereas HIC15 and HARM showed weaker correlation with MPS, particularly in rugby. XGB showed the lowest cross-sport variability (CV=0.034) and uncertainty (median NCIS=0.056). HARM, DAMAGE and HIC15 showed the greatest sport dependence (CV=0.575-0.588) and uncertainty (median NCIS=0.331-0.791). Conclusions: XGB, BrIC, and UBrIC demonstrated the strongest associations with brain tissue strain and the greatest consistency across sports. This study provides a biomechanically informed framework for selecting suitable metrics for sports HIA protocols.

Cells adjust their internal circuits in response to changes in their environment. Hence, exposing cells to changing conditions provides a way to probe the intrinsic dynamics of cellular internal circuits. Metabolic networks are examples of such circuits since metabolic fluxes dynamically adjust when environmental conditions are transiently altered. Most existing theoretical frameworks focus on cellular metabolic steady states and do not consider the dynamics of changes in metabolic fluxes. In this work, we applied transfer function analysis from control theory to analyze the changes of NADH oxidative fluxes in the mitochondria and cytoplasm in mouse oocytes in response to dynamical perturbations of oxygen depletion and recovery. We observed an overshoot of NADH oxidative flux in the cytoplasm upon oxygen recovery which is absent in the mitochondrial NADH oxidative flux. Metabolic perturbation experiments and transfer function analysis indicate that this cytoplasmic NADH overshoot results from the coupling of the mitochondrial and cytoplasmic NADH cycles. The degree of overshoot is determined by competing timescales associated with the exchange rates of lactate and pyruvate with the media and their interconversion rates catalyzed by lactate dehydrogenase. Applying control theory to the data enables the inference of the exchange and conversion rates of pyruvate and lactate, allowing predictions of the contribution of lactate to mitochondrial respiration. Our work indicates that the oocytes maintain a homeostatic respiration rate across nutrient conditions by modulating the contribution of lactate to mitochondrial respiration.

Background: Acute pulmonary embolism (PE) is a leading cause of morbidity and mortality with persistent difficulties in choosing interventions and predicting outcomes for patients defined clinically as intermediate risk. Computational fluid dynamics (CFD) tools have been used to understand the hemodynamic environment and plan interventions in the pulmonary arteries across a variety of disease conditions. Several biomechanical metrics have been used to evaluate risk in narrowed vessels, including hemodynamic resistance, power dissipation, and fractional flow reserve (FFR). In this study, we evaluate differences in these CFD-derived biomarkers between healthy controls (HC) and intermediate risk, acute PE patients. Additionally, we examine the response of patient hemodynamics to mechanical thrombectomy and compare values of these biomarkers across post-intervention pressure status. Methods: A CFD framework was developed to simulate patient-specific hemodynamics within the pulmonary vasculature identifiable from clinical imaging. The pipeline involved reconstructing three-dimensional (3D) structures of the pulmonary arteries and modeling blood flow with the finite element method. Patient-specific boundary conditions were derived from matching pre-intervention inlet mPAP to the patient's measured value given their measured CO as steady inflow. Converged simulations allowed for precise quantification of primary hemodynamic characteristics (flow and pressure) as well as secondary flow phenomena, primarily wall shear stress (WSS) and simulated pressure metrics such as fractional flow reserve (FFR). Results: Our simulations revealed significant elevations in resistance, power dissipation, and the number of vessels with low FFR in those patients with acute PE (n=6) compared to HC (n=3). Occlusions of hemodynamic significance were generally found in segmental pulmonary arteries. For patients with normalized pulmonary pressures post-thrombectomy (n=3), we found significantly higher proximal power dissipation and counts of low FFR vessels in comparison to those with elevated pressures after intervention (n=3). Distal resistance, which was derived from the portion of resistance attributed to the outflow boundary conditions, was significantly higher in patients with elevated pressures post-intervention. Across all PE patients, FFR count was significantly correlated with post-thrombectomy pulmonary pressure and cardiac index. Discussion: CFD-derived biomarkers offer a promising tool for understanding disease severity in acute PE. Differences between HCs and acute PE patients reveal expected increases in metrics associated with proximal disease burden. Yet, in examining acute PE patients with varying post-intervention hemodynamics, we found that these metrics of proximal disease burden could also be useful to predict the efficacy of mechanical thrombectomy. Those patients with normalized pressures had higher values for proximal disease metrics and lower values for distal disease metrics than those with continued elevations in pressure. This suggests that accessibility of hemodynamically-significant emboli to thrombectomy may be useful as a predictor for outcomes.

We present three quantitative methods: 1) estimation of exoskeleton mechanical power and energy ratio from published data, 2) a systematic review of the exoskeleton literature on reported energy ratios, and 3) timing correction analysis of the replication experiment, to address concerns raised by Collins et al. (2026) about Luo et al. (2024). Together, these analyses support the reported metabolic reductions and the validity of exoskeleton control via learning in simulation. The critique rests on an unsupported premise: that exoskeleton energy ratios above 4 are physiologically implausible. This premise of Collins et al. (2026) is not supported by the cited evidence, and the error originates in their own cited source. Sawicki and Ferris (2009), the paper they invoke as authority for the limit of 4, state explicitly that "reported values of the muscular efficiency range from 0.10 to 0.34, with many sources assuming an average of [~]0.25." The value of 4 corresponds to this average, it is not a physiological ceiling. Treating an average as a physiological upper limit is a fundamental error. The published exoskeleton literature further contradicts the claim, including work by the authors of the critique themselves (Collins et al., 2015: 4.3; Young et al., 2017: 5.0) and independent work (Malcolm et al., 2013: 4.8; Seo et al., 2017: 6.7). In contrast, our walking energy ratio is 2.4, calculated directly from Fig. 4 of our paper. Our device delivers higher peak torque (14.1 Nm vs. 10.9 Nm, Lim et al., 2019) and achieves a slightly larger metabolic reduction (24.3% vs. 21%). Independent groups have since demonstrated meaningful metabolic reductions using learning-in-simulation frameworks, including Barati et al. (2026, 15.2% mean and 22.5% maximum) and Zhou et al. (2025, [~]20% during running). The claim of Collins et al. (2026) that this problem "remains unsolved" is directly contradicted by these independent results. The experiment in the critique is not a valid replication of our method. Our controller is a neural network with [~]10,000 parameters learned through deep reinforcement learning in musculoskeletal simulation; the critique instead applies a pre-programmed fixed torque curve with no learnable parameters. Beyond this, the replication contains three methodological errors: 1) a heel-strike timing assumption producing offsets up to 30% of the gait cycle; 2) an averaged torque profile that discards subject-specific control; and 3) a device [~]50% heavier than ours (4.8 kg vs. 3.2 kg) without measuring the metabolic penalty of the added weight. The critique also misreports Samsung data, with reported values approximately double those in the original publication, errors that directly underpin their physiological limit argument.

PurposeAnkle stiffness can be altered by normal aging, bone and joint pathology, and treatments such as orthoses or surgical joint fusion. The effects of ankle stiffness on gait are not yet well understood but may be crucial for understanding how these pathologies and treatments influence body mechanics. The objective of this work was to investigate how isolated changes in stiffness applied in parallel with the ankle impact lower-limb kinematics, kinetics, joint work, and muscle activation during walking in individuals without lower-limb pathology. MethodsNine young adults without lower-limb pathology wore an adjustable-stiffness ankle exoskeleton and walked at 31 different conditions of ankle spring stiffness, neutral angle, and treadmill incline. We recorded motion capture data, ground reaction forces, and muscle activation, and analyzed the resultant data for trends as a function of ankle stiffness. ResultsExoskeleton-side ankle range of motion decreased and asymmetry increased across all joints as ankle stiffness increased, primarily due to decreased plantarflexion at toe-off. The 30 Nm/rad spring stiffness condition led to a minimum in mean exoskeleton-side muscle activation and hip joint work, but increased kinematic asymmetry. ConclusionOur results suggest that there may exist a range of stiffnesses at the lower end of typically-studied values that can reduce muscle activation and joint work during walking, though at the cost of kinematic symmetry. These findings provide a deeper understanding of how ankle stiffness influences gait mechanics, with potential applications in wearable devices, clinical rehabilitation, and assistive technology.

Cetaceans reduce the net cost of sustained swimming through intermittent locomotion, alternating active fluking with unpowered gliding. The energy balance of this strategy is central to understanding survival rates, population sustainability, and the effects of anthropogenic and environmental pressures. While active-phase energetics have been characterized extensively, the glide phase remains largely unexplored. Here we derive the optimal glide duration (Topt) and the maximum glide duration beyond which energy savings vanish (Tzero) for three odontocetes spanning a 20-fold range in body mass, using high-fidelity CAD models and wall-modeled large eddy simulations. We show analytically that speed retention at Topt and mass-specific peak energy savings are both fully determined by the active-to-passive drag ratio, propulsive efficiency, and swimming speed, independently of body morphometry and drag coefficient, and are therefore invariant across species at any given speed. These passive-phase optima extend the known size-independent active-phase invariants to the glide phase, towards a scale-independent energetic framework for burst-and-glide locomotion in small cetaceans.

Energy expenditure (EE) during running depends on the interplay between active muscle work and elastic energy storage and return, yet the relative contribution of mechanical power to EE remains debated. Quantifying the relative contributions of segment-level mechanical power can provide a way to address this debate. In this study, we aimed to quantify how segment-level mechanical power contributes to EE during running and to demonstrate that these mechanistic insights support wearable-based EE estimation. Joint dynamics and respiratory gas-based EE were collected from healthy young adults running at multiple speeds. Scale factors were derived to quantitatively link efficiency-weighted segment power to measured EE. The stance leg consistently showed the strongest correlation with EE, and this dominance was preserved across speeds. Including swing-leg hip power further improved accuracy. Scale factors were approximately 0.45, suggesting that active muscle work and elastic energy return contribute comparably to the mechanical power associated with EE. Using a lightweight machine learning model, stance-leg and swing-leg hip joint power were reconstructed from a single sacral IMU, enabling accurate EE prediction. These findings demonstrate that lower-limb mechanical power is a robust predictor of running EE, supporting both the extensibility of biomechanically-informed frameworks and wearable-based EE monitoring.

ObjectiveOur prior studies have demonstrated that lateral spinal cord stimulation can evoke somatosensory percepts in the missing foot in individuals with a lower-limb amputation. However, subjects reported concurrent sensations in their residual limb. In this study, we evaluate the hypothesis that using high-density paddle electrodes with smaller contact sizes, and multipolar stimulation configurations could evoke more focal sensations in the foot over a wide range of stimulation amplitudes. ApproachWe used a combination of electrophysiology and computational modelling methods to investigate the selective activation of distal nerve branches in response to lateral spinal cord stimulation in cats. In six acute feline experiments, we performed an L3-S1 laminectomy and placed custom 32-electrode paddles laterally over the dura of the spinal cord. We recorded antidromic action potentials in the distal branches of the sciatic and femoral nerve trunks in response to stimulation using three contact diameters (150, 500 and 1000 {micro}m) and two stimulation configurations - monopolar and bipolar stimulation. We replicated the neural recruitment patterns from those experiments in a computational model of the feline lumbar spinal cord. We then used the model to examine neural recruitment with 1.8 mm and 2.5 mm contacts, as well as a tripolar guarded-cathode configuration. Main resultsIn the electrophysiology experiments, the 500 {micro}m-diameter electrodes achieved the most selective nerve activation (68%) compared to 62% for both 150 and 1000 {micro}m-diameter electrodes. The minimum amplitudes for recruiting nerve branches (i.e., threshold) as well as the dynamic ranges were largely similar for the different contact diameters (median: 35 {micro}A) and stimulation configurations (30 {micro}A for bipolar stimulation; 35 {micro}A for monopolar stimulation). The computational model reproduced the finding that selectivity did not differ significantly among the three contact sizes tested in cat experiments, though it revealed that increasing contact diameter above 1000 {micro}m raised the minimum amplitude required for selective activation and reduced spinal root selectivity. Across both approaches, we consistently recruited large-diameter afferents that are critical for somatosensory applications of spinal cord stimulation. SignificanceOur results indicate that, relative to clinical electrodes, reducing the contact diameter of stimulation electrodes can evoke focal sensations, but further reductions below 1000 {micro}m may fail to improve selectivity. This study highlights potential constraints with achieving focal selectivity that are not dependent on the design of the electrodes.

The endothelial cells (ECs) that line blood vessels continuously sense and respond to the physical forces exerted by blood flow. In vivo, pulsatile arterial flow interacts with vessel curvature, branching and other anatomical features to generate complex local hemodynamic environments that dictate the magnitude, direction, pulsatility, and oscillatory nature of wall shear stress experienced by ECs. Currently, accessible and reproducible in vitro models of complex pulsatile flow that recapitulate in vivo vascular anatomy remain limited. Here, we combine a novel rotational-flow endothelial culture platform with detailed computational fluid dynamics (CFD) modeling to characterize four well geometries designed to generate distinct hemodynamic environments. CFD analyses demonstrate that these geometries intrinsically generate pulsatile flow and produce reproducible spatially distinct regions of wall shear stress magnitude, pulsatility, and oscillatory shear within a single culture well. Endothelial alignment mapping and functional assays reveal region-specific cellular responses to the predicted local flow conditions that closely corresponded to the predicted local hemodynamic environment, linking complex flow patterns to endothelial adaptation. The technical advancements of our modeling efforts should support a faster, cheaper, simpler, and--importantly--validated framework for future investigation into EC mechanobiology under complex flow conditions. HIGHLIGHTSO_LISimple engineered well geometries generate distinct hemodynamic microenvironments, mimicking in vivo vascular structures, using a conventional orbital shaker. C_LIO_LIComputational fluid dynamics (CFD) reveals spatially distinct patterns of wall shear stress, pulsatility, and oscillatory shear applied to ECs within individual culture wells. C_LIO_LIHigh average wall shear stress and elevated oscillatory shear index induces a unique perpendicular alignment of ECs to the dominant flow vector. C_LI

The human masseter muscle is one of the primary muscles responsible for mastication and mandibular movement; however, its intrinsic mechanical properties remain insufficiently characterized. In this experimental study, the nonlinear, viscoelastic, and history-dependent behaviour of the human masseter muscle was investigated using ex vivo uniaxial cyclic tensile tests. The masseter muscle samples prepared from fresh and formalin-preserved cadavers were tested under two loading protocols: a continuous stretch protocol with increasing stretch levels and a constant stretch protocol with repeated loading to a fixed maximum stretch. Tests were conducted at two strain rates, and their influences on the mechanical behaviour of the tissue were examined. The effect of formalin preservation was also investigated. The results showed that the stiffness of the tissue increases for formalin-preserved samples. Under cyclic loading, the features including energy dissipation, stress-softening, residual deformation, and cyclic conditioning progressively changed during the initial loading cycles and reached stabilization during the final cycle. These findings provide experimental evidence that the human masseter muscle exhibits nonlinear, viscoelastic, and history-dependent mechanical behaviour under cyclic tensile loading. The experimental data obtained in this study may be used for biomechanical modelling of the human masticatory system and the development of constitutive models for cranio-maxillofacial surgical simulation, prosthetic design, and facial soft-tissue biomechanics. Statement of significanceThe masseter muscle is one of the primary muscles of mastication. To address the current gap in craniofacial biomechanics that has largely focused on the mechanical characterization of the masseter muscle based on imaging techniques or monotonic loading, this study quantifies the nonlinear and viscoelastic mechanical response of masseter tissue under cyclic continuous and constant stretch loading, including strain-rate and preservation effects. The results show that the mechanical behaviour of the masseter muscle, including stiffness, hysteresis, stress-softening, and residual strain behaviour, is strongly influenced by strain-rate and formalin preservation. The experimental results provide mechanical data for constitutive modelling of the masticatory system with applications in cranio-maxillofacial surgical simulation, prosthetic design, and facial soft tissue modelling.

Accurate modeling of laryngeal motor control is key to understanding typical and disordered voice production. However, traditional biomechanical plant models based on ordinary differential equations (ODEs) often involve high computational costs and numerical instabilities, limiting their use in real-time closed-loop control frameworks. This study evaluates feature-driven machine learning (ML) regressors, specifically Random Forest (RF), Multilayer Perceptron Neural Networks (NN), and Polynomial Regression (PR), as surrogate forward models mapping laryngeal motor inputs to fundamental frequency and sound pressure level. Training data were generated with two biomechanical vocal fold models: the extended body-cover and the triangular body-cover. Results demonstrate that ML surrogates reduce execution times from seconds to milliseconds (e.g., 2 ms for PR), enabling stable real-time tracking via inverse Jacobian control. While RF provides the highest accuracy, NN and PR offer smoother control signals and smaller memory footprints. A practical performance threshold was identified near N = 1,000 training samples, below which accuracy degraded substantially when models were trained from scratch. These findings support ML surrogates as efficient and adaptable alternatives to direct numerical simulation, providing a foundation for future subject-specific modeling through transfer learning in data-limited clinical scenarios.

The development and physiology of neurons rely on their microtubule organization, which is characterized by plus-end-out oriented microtubules in the axon and a mix of plus-end-out and plus-end-in oriented microtubules in dendrites. This orientational pattern is established early in neuronal development and is tightly linked to axon-dendrite differentiation. Even though multiple potentially relevant mechanisms have been proposed, fundamental questions remain: How does the microtubule organization in neurons emerge, and how does a neuron develop a single axon and multiple dendrites? Here, we address these questions at two distinct, complementary levels: at a higher level by proposing a conceptual framework, in which we classify mechanisms into three categories based on how they contribute to the microtubule organization: orientational bias, parallel amplification, and polarization; at a lower level we build a biophysical model that incorporates multiple mechanisms of microtubule dynamics in a neuron, from which, using analytical calculations and simulations, we derive insights into the emergence of microtubule organization in developing neurons. We show that geometrical effects alone can confer a bias in microtubule orientation. Parallel amplification then enhances the resulting polarity. Coupling multiple neurites to a common cell body that serves as a shared reservoir of resources allows for a polarization mechanism that ensures that the microtubule organization of one neurite becomes axonal while all others are dendritic. This framework unifies diverse molecular observations and yields experimentally testable predictions about microtubule self-organization in early neuronal development. Author summaryNeurons communicate through long protrusions called neurites, which are of two types: dendrites, which receive signals, and axons, which send signals. Their development relies primarily on microtubules, which are polar filaments with two distinct ends, known as the plus and minus ends. Microtubules self-organize into functional architectures that are significantly different between axons and dendrites. In axons, all microtubules point their plus end away from the cell body, whereas in dendrites, they point either towards the cell body or have mixed orientations depending on the species. This orientation guides intracellular transport by motors and is closely linked to whether a neurite develops into an axon or a dendrite. Despite decades of research identifying individual mechanisms, the bigger picture behind the emergence of microtubule orientation in neurons remains unclear. Here, we construct a conceptual framework and a biophysical model to identify the principles underlying the emergence of microtubule orientation in developing neurons. Our conceptual framework provides a high-level perspective on how individual mechanisms influence microtubule organization in neurites. In our concrete biophysical model, we study a selection of mechanisms to gain specific, quantitative insight into the organizational process. We propose a minimal model of a neuron that exhibits neuronal polarization, giving rise to a single axon-like neurite and multiple dendrite-like ones, consistent with experimental observations. This in silico neuron helps to explain how neurons break symmetry during development and provides a systematic way to generate and test new hypotheses about neuronal polarity.